from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C+2")


id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 6.5,
    A   = -4.024e-2,
    B   = 2.914,
    C   = -4.344,
    D   = 2.121,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s2p3P",
    X1   = 2.0,
    P   = 5.434e-1,
    Q   = 3.453e-1)


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 12.7,
    A   = 6.261e-1,
    B   = 3.044,
    C   = 5.384e-1,
    D   = 0.0,
    E   = 4.364,
    F   = 5.769e-2,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s2p1P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 17.2,
    A   = -5.732e-3,
    B   = 8.467e-2,
    C   = -1.676e-1,
    D   = 1.307e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2p23P",
    X1   = 2.0,
    P   = 5.389e-3,
    Q   = 7.167e-3)


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 18.2,
    A   = 2.915e-1,
    B   = 1.655e-1,
    C   = -3.445e-1,
    D   = 3.720e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2p21D")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 22.9,
    A   = 2.092e-2,
    B   = 2.289e-1,
    C   = -4.268e-1,
    D   = 2.137e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2p21S")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 29.5,
    A   = 3.910e-3,
    B   = 2.740e-1,
    C   = 1.120,
    D   = 1.110,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3s3S")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 30.6,
    A   = 4.260e-1,
    B   = 4.350e-1,
    C   = -2.280,
    D   = 1.760,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3s1S")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 32.1,
    A   = -1.140e-1,
    B   = -5.230e-1,
    C   = 8.110e-1,
    D   = 0.0,
    E   = 3.730e-1,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3p1P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 32.2,
    A   = -2.190e-2,
    B   = 3.660e-1,
    C   = -1.030,
    D   = 8.240e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3p3P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 33.5,
    A   = 4.190e-2,
    B   = 6.490e-2,
    C   = 1.100e-2,
    D   = 1.050e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3d3D")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 34.3,
    A   = 7.330e-1,
    B   = 5.400e-2,
    C   = -2.110,
    D   = 1.490,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+2_2s21S-2s3d1D")

cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()